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ABSTRACT 

Recognition  of  class  and  aspect  angle  of 
a  rigid  three-dimensional  object  in  free  space 
is  facilitated  by  using  N-dimensional  chain 
codes  in  feature  space.  In  both  radar  and  im¬ 
age  recognition  systems  features  have  previ¬ 
ously  been  described  that  are  invariant  to  ob¬ 
ject  rotation  about  the  observer's  line  of 
sight  to  the  object  center  of  gravity.  Howev¬ 
er,  rotations  about  the  remaining  two  orthogo¬ 
nal  axes  produce  changes  of  the  features  and  a 
full  description  of  the  object  becomes  a 
toroidal  surface  in  N-dimensional  feature 
space.  N-dimensional  chain  codes  are  intro¬ 
duced  in  this  paper  to  give  a  line  approxima¬ 
tion  of  such  a  surface  that  is  easily  used  to 
identify  an  unknown  object  and  to  estimate  its 
aspect  angle. 

1-  INTRODUCTION 

N-dimensional  chain  codes  are  a  direct 
extension  of  the  familiar  and  extremely  useful 
concept  of  two-dimensional  chain  codes,  which 
are  used  to  make  piecewise-linear  approxima¬ 
tions  of  continuous  curves  with  a  small  set  of 
directed  line  segments.  In  this  paper 
straight-forward  procedures  are  given  for  gen¬ 
erating  a  chain  code  in  N-dimensions,  either 
from  a  continuous  curve  or  from  sample  points 
of  the  curve.  A  harmonic  description  of  a 
closed  chain  code  is  derived  as  well  as  a 


bound  on  the  error  introduced  by  truncating 
the  harmonic  content.  Comparisons  are  made  of 
the  information  required  to  encode  a  curve 
with  a  chain  code  versus  conventional  full- 
coordinate  description  of  approximation  points 
on  the  curve.  N-dimensional  codes  can  be  used 
for  efficient  storage  of  toroidal  classifica¬ 
tions  in  an  N-dimensional  feature  space  and  to 
facilitate  comparison  of  known  classifications 
with  the  features  of  an  unknown  object.  An 
application  of  N-dimensional  codes  is  shown  to 
the  successful  recognition  of  class  and  aspect 
angle  of  aircraft. 

2.  CHAIN  ENCODING  0£  2-DIMENSIONAL  ARBITRARY 
6E0METRIC  CURVES 

A  method  of  angular  quantization  and  di¬ 
gital  encoding  of  arbitrary  2-dimensional 
geometric  curves  has  been  proposed  by  Freeman 
and  is  extended  here  for  N-dimensional  curves. 
The  method  for  2-dimensional  curves  consists 
of  superimposing  a  square  grid  on  the  curve, 
and  approximating  the  curve  with  those  grid- 
intersection  points,  in  the  order  of  a  trace 
progressing  along  the  curve,  which  lie  closest 
to  the  curve.  An  example  of  a  curve  approxi¬ 
mated  in  this  manner  is  shown  in  Figure  1.  A 
trace  starting  from  point  X  is  moved  along  the 
curve  until  a  grid  line  is  intersected  at 
point  B.  The  distances  AB  and  CB  are  com¬ 
pared,  and  since  AB  is  less  than  CB,  the  grid 
intersection  point  A  is  chosen  as  an  approxi- 


Fig.  1.  Approximation  of  a  curve  by  closest  grid  intersection  points. 
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nation  point.  The  trace  then  proceeds  to  the 
next  grid  intersection  at  point  0.  The  dis¬ 
tance  ED  is  less  than  CD;  therefore  point  E  is 
chosen  as  an  approximation  point.  The  trace 
intersects  a  grid  line  again  at  point  F,  and 
EF  is  less  than  GF,  but  point  E  has  already 
been  designated  as  an  approximation  point  and 
cannot  be  used  twice  in  succession.  The  ap¬ 
proximation  points  H,  I,  J  are  found  in  a 
similar  manner  as  the  trace  progresses  along 
the  curve  and  intersects  grid  lines  at  points 
L,  M,  N,  respectively.  If  A8  =  C8,  for  exam¬ 
ple,  a  consistent  choice  of  approximation 
points  is  obtained  by  using  the  candidate 
point  closest  to  the  origin  of  the  coordinate 
axes.  Since  the  curve  is  assumed  to  be  con¬ 
tinuous,  it  can  progress  from  an  approximation 
point  to  only  one  of  eight  possible  neighbor¬ 
ing  points  on  the  grid.  Each  possible  pro¬ 
gression  between  a  pair  of  approximating 
points  is  called  a  link  or  vector,  and  the 
curve  is  described  by  the  succession  of  these 
vectors.  For  example,  the  vector  chain  for 
the  curve  in  Figure  1  is  shown  by  the  sequence 
of  four  straight  dotted  lines  between  points 
A,  E,  H,  I,  J. 

Each  vector  encodes  an  incremental  change 
in  position,  length,  and  directional  angle  of 
the  trace  as  it  progresses  along  the  curve. 
It  is  proposed  here  that  the  chain  code  con¬ 
sist  of  the  sequential  listing  of  the  incre¬ 
mental  changes  of  position  (Ax, Ay)  for  each 
vector;  e.g.,  the  link  chain  V.,  in  Figure  1  is 

as  follows: 

V,  =  <1, 1)0, OCT, 0)  (1,-1 ) 

The  incremental  change  in  length  due  to  a  par¬ 
ticular  vector  (Ax, Ay)  is  At  =  (Ax^+Ay^)^^ 
and  the  directional  cosines  of  the  vector  for 
the  x  and  y  dimensions  are,  Ax/At  and  Ay/At, 
respectively.  Note  that  Ax  and  Ay  can  take  on 
only  the  values  -1,  0,  and  1. 

O  CHAIN  ENCODING  N-DIMSNS IONAL  ARBITRARY 
GEOMETRIC  CURVES 

A  continuous  arbitrary  N-dimensional 
curve  can  be  chain  encoded  by  mapping  to  in¬ 
tersection  points  of  an  N-dimensional  grid. 
Approximation  points  are  generated  in  a 
fashion  completely  analogous  to  the  2-dimen¬ 
sional  case  when  the  curves  crosses  a  quanti¬ 
zation  interval.  The  crossing  point  is  mapped 
to  the  nearest  grid  intersection  and  a  vector 
V ,  is  generated,  consisting  of  the  incremental 

variation  in  each  dimension  between  the  ap¬ 
proximation  points  P,_,  and  P..  The  vector 

(V)  is  expressed  as 

N* 

V  s  T  v.  (1) 

i*1 


where  v,  has  N  components  (X,,,  X^,  ...,  X-N> 

corresponding  to  the  number  of  dimensions  in 
the  space.  A  3-dimensional  example  of  vector 
labeling  and  encoding  is  shown  in  Figure  2. 
Note  that  since  the  maximum  variation  between 
approximation  points  is  one  grid  unit  in  each 
dimension  a  vector  component  X.^  can  only  take 

on  the  values  of  -1,  0,  and  +1.  A  vector 
train  in  N  dimensions  may  be  projected  into  a 
lower  L-dimensional  space  by  omitting  the  com¬ 
ponents  of  the  original  vector  in  the  N-L  di¬ 
mensions  not  of  interest  and  eliminating  any 
resulting  null  vectors.  In  Figure  2,  the  pro¬ 
jection  of  the  vector  train 


V123  =  11 ,0,-D  (0,1,0)  (0,1, OHO,  1,-1)  (1,0,0) 

(1,-1 ,0)  (0,-1 ,1)  (1,0,0) 

(1,-1)  (0,-1 )  (1 ,0)  (1 ,0)  (0,1 )  (1 ,0) , 

where  two  null  vectors  have  been  removed. 

The  length  At.  of  the  vector  v,  is 

At.  =  (£  |X..|r'*  (2) 

1  £i 

and  the  length  tq  of  the  first  q  links  in  the 
chain  is 

tq  »  |  At, 


H 


Fig.  2.  A  3-D  vector  chain  approximation  of  a 
curve  and  its  projection  on  the  x,-Xj  plane. 
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For  a  closed  vector  consisting  of  K  vectors 
the  period  T  of  the  chain  is 

T  =  lK  U) 
4^  INFORMATION  CONTENT  0F_  THE  CHAIN  CODE 

For  N  dimensions  there  are  3  possible 
vector  types,  including  the  null  vector,  re¬ 
quiring  a  minimum  of  N  log2  3  binary  bits  of 

computer-storage  per  vector.  The  correspond¬ 
ing  bit  storage  requirement  for  a  single  ap¬ 
proximation  point  is  N  log^  S,  where  S  is  the 

value  of  the  maximum  coordinate  in  the  space. 
The  ratio  of  total  bit  storage  requirements 
for  recording  the  full  coordinate  information 
for  K  approximation  points  along  a  continuous 
curve  versus  vector  encoding  for  the  same  K 
approximation  points  is 


found  because  x.(t>  is  piecewise  li near  and 

continuous  for  all  time.  The  derivation  of 
the  coefficients  here  involves  the  time 


derivative  x.(t),  which  consists  of  the  se¬ 
quence  of  piecewise  constant  derivatives 

Ax  ./At  for  values  of  p  in  the  range  of 
PI  P 

1  <  p  K.  The  time  derivative  is  also 
periodTc  with  period  T  and  by  using  this  ob¬ 
servation  it  is  easily  concluded  by  Kuhl  and 

Giardina^  that  the  Fourier  coefficients 

a.  and  b.  are  as  follows: 
in  )n 


A. 


i°  =  T  2itP 


p  VI 5  +  WVi5 


(7) 


where. 


K  N  log2  S 

CK-1 T  N_ log2  3  +  N  log2  S 


(5) 


P-1 

=  £ 

i=1 


Axpj  "  nf"Ati 


For  Large  K,  the  ratio  approaches  log^  S  and 

vor  values  of  S  >  3  the  vector  encoding  pro¬ 
vides  significant  economy  in  storage  require¬ 
ments. 

The  minimum  bit  requirement  for  a  vector 
can  be  approached  only  by  using  special  coding 
schemes.  The  particular  cases  for  N  =  2  and 

N  =  3  have  been  treated  by  Freeman1  and  Rut- 
2 

tenburg  .  For  N  ^  4,  since  each  component  can 
take  on  3  values,  and  because  of  the  binary 
nature  of  digital  computers,  2N  bits  per  vec¬ 
tor  are  often  used  and  the  storage  ratio  will 
approach  .5  log2  S  =  .793  log^  S  for  large  N. 

For  this  more  general  case,  storage  economy  is 
realized  when  S  >  4. 

5 .  HARMONIC  DESCRIPTION  0f_  A  CLOSED 
N-DIMENSIONAL  CHAIN  CODE 

The  Fourier  Series  expansion  of  the  pro¬ 
jection  of  a  closed  N-dimensional  chain  code 
on  the  Xj  axis  is  defined  as 

V°  *  V£  ain cos  ^bjn sin  ^r1  <6> 

where, 

*jo  *  t/o  *j(t>  dt 

•,n-f/{  »j‘«  cos  ^11  dt 
bjn  ’  7/5  *j(t>  sin^dt 


and. 

T 

K 

&xpj 

2nxt 

cn-  0 

ajn 

2n2*2 

2- 

P=1 

«p 

cos  T 

T 

K 

v 

A>tpj 

2nxt 

e  sn  * 

bjn 

2n2*2 

2- 

P=1 

% 

sin  j 

2n»t 


P-1 


(8> 


2n»t 


-  sm 


P-1 


(9) 

It  is  useful  to  be  able  to  specify  the 
number  of  harmonics  M  required  such  that  a 
truncated  Fourier  approximation  to  the  projec¬ 
tion  Xj  of  a  closed  vector  chain  have  a  max¬ 
imum  absolute  error  no  greater  than  e^.  Let 


JN 


=  Aj0*£.  ai 


n=1 


2nwt  .  .  2nwt 

jn  cos  "T  +  bjn  Sln  —  (10) 


It  is  shown  by  Giardina  and  Kuhl 
bounded  by  the  expression 


that  t . 
J 


2»  N 


•  V  CXjCt)) 


(11) 


For  a  closed  vector  chain 


I  <;,»» .  h  |ipj-v,j|  (iz> 


The  Fourier  coefficients  corresponding  to  the 

nth  harmonic  and  b.  (b,„  s  0)  are  easily 
jn  jn  JO 


where. 


3 
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The  error  of  the  Fourier  approximation  to  the 
closed  chain  code  truncated  after  M  harmonics 
is  then  bounded  by  the  maximum  Cj. 

Harmonic  descriptions  derived  from  chain 
codes  for  2-dimensional  closed  figures  can  be 
normalized  in  straight-forward  fashion  to  be 
independent  of  size,  translation,  rotation  and 
dilation  and  the  particular  set  of  descriptors 

of  Wallace  and  Wintz*  are  used  in  a  later  sec¬ 
tion  of  this  paper  to  obtain  experimental 
results  for  aircraft  recognition.  Another 
such  set  of  descriptors  was  developed  by  Kuhl 

and  Giardina^.  It  is  envisioned  that  similar 
harmonic  descriptors  could  be  developed  for 
higher  dimensional  closed  curves  by  using  the 
elliptic  properties  of  the  Fourier  coeffi¬ 
cients  described  by  Kuhl.* 

6 .  FOURIER  DESCRIPTION  0£  THREE-DIMENSIONAL 
OBJECTS 

This  section  concerns  the  identification 
of  three-dimensional  rigid  objects  when  viewed 
by  a  two-dimensional  video  imaging  system  from 
an  arbitrary  aspect  angle. 

If  a  rigid  object  in  free-space  is  imaged 
using  visible  or  somewhat  longer  wavelengths 
the  resulting  image  can  be  digitized,  the 
shape  (usually  silhouette  or  contour)  extract¬ 
ed,  and  the  object  recognized.  Two  common 
features  used  are  two-dimensional  measurements 

of  the  silhouette  or  contour7-1*^  and  Fourier 

descriptors  of  the  contour. ^ 

The  method  used  in  the  experiments 
described  in  the  later  section  is  Fourier 
descriptors  as  discussed  in  Sec.  5  of  this  pa¬ 
per  fully  desrribed  in  References  5  and  15. 
These  features  can  be  normalized  to  remove  ef¬ 
fects  of  object  scale,  translation,  and  rota¬ 
tion  about  the  imaging  axis  (9^  in  Fig.  3). 

However  two  degrees  of  freedom  still  remain 
for  movement  of  the  three-dimensional  object 
which  correspond  to  change  of  viewing  aspect 

angles  (0  and  9  in  Fig.  3).  For  an  arbi- 
x  y 

trary  three-dimensional  object,  viewing  aspect 
angle  changes  will  create  different  two- 
dimensional  representations  at  the  sensor  and 
each  view  must  be  stored  as  a  separate  library 
entry  for  classification  purposes.  As  an  ex¬ 
ample,  Fig.  4  shows  77  views  of  an  airplane 
representing  various  aspect  angles.  Those  li¬ 
braries  may  be  stored  using  the  N-dimensional 
chain  code  concepts  described  earlier. 

The  normalized  Fourier  descriptor  ob¬ 
tained  consists  of  coefficients  representing 
harmonics  of  the  contour.  Usually  at  least  16 
harmonics  are  required  (64  reat  numbers)  to 


Fig.  3.  Coordinate  system  for  imaging  three- 
dimensional  object. 
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Fig.  4.  2-D  views  resulting  from  the  aspect 
angles  used  for  storing  the  library  of  an  F105 
airplane.  Horizontally  are  8^  angles  ranging 

from  0  to  x/2  radians,  and  vertically  are  8 
angles  ranging  from  -  x/2  to  +  »/2  radians. 
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adequately  represent  a  complex  contour.  How¬ 
ever  for  a  given  class  of  problems  such  as 
airplane  recognition,  considerable  correlation 
exists  among  these  coefficients.  Thus  an  ap¬ 
propriate  linear  transform  (based  on  eigenvec¬ 
tors  of  the  autocorrelation  matrix  of  all  pos¬ 
sible  NFDs)  is  used  to  reduce  the  dimensional¬ 
ity  of  the  data.  Approximately  12  real  num¬ 
bers  must  be  retained  to  recognize  various 
airplanes  from  different  aspect  angles. 

Thus  if  a  rigid  object  in  free  space  is 
imaged  and  the  resulting  contour  converted  to 
a  normalized  (N-dimensional)  Fourier  Descrip¬ 
tor,  the  effects  of  rotations  are  normal¬ 
ized  out,  but  changes  in  8  and  8  give  other 

x  y 

aspect  views  whirh  give  different  N- 
dimensional  vectors. 

7_.  DIFFERENTIAL  LIBRARY  REPRESENTATION 

In  order  to  reduce  storage  requirements 
and  to  speed  the  classification  process,  a  new 
method  of  storing  the  library  information  has 
been  investigated.  He  are  considering  the  li¬ 
brary  of  each  object  as  a  set  of  sampled 
points  from  an  N-dimensional  surface  (N 
represents  the  number  of  features  retained,  12 
in  our  examples).  This  surface  represents  all 
possible  N-dimensional  vectors  obtained  as  the 
aspect  angle  is  changed.  It  has  the  topology 
of  a  toroid.  If  we  code  the  N-dimensional 
vectors  independently,  we  require  k'N  bytes 
for  each  vector  where  k  is  the  number  of  bytes 
used  to  code  each  feature  value. 

Instead,  we  propose  to  arrange  the  sample 
vectors  in  some  order  and  code  the  N- 

fimensional  difference  between  successive  sam- 
;s.  If  the  difference  in  each  dimension  is 
only  allowed  to  have  values  of  -S,  0,  or  +S  we 
have  a  form  of  N-dimensional  delta  modulation 
with  step  size  S.  The  library  can  be  stored 
using  the  N-dimensional  code  described  in  Sec¬ 
tion  III. 

Three  potential  advantages  are  possible 
using  such  a  scheme: 

(1)  if  the  vectors  are  arranged  in  an  order  so 
that  successive  vectors  are  similar 
(correlated)  the  differential  coding 
scheme  may  allow  representation  of  the  N- 
dimensional  surface  using  less  storage; 

(2)  if  the  distance  measure  can  be  incremen¬ 
tally  calculated  using  the  differential 
vector  information.  As  the  vectors  are 
sequentially  compared,  the  time  required 
to  search  the  library  may  be  reduced; 

(J)  when  the  step  size  S  is  small,  the  inter¬ 
mediate  vectors  generated  between  library 
samples  provides  a  form  of  interpolation 
which  is  much  faster  than  typical  interpo¬ 
lation  methods. 

The  potential  disadvantages  of  such  a 
method  are: 

(1)  the  quantization  error  introduced  by  a 


fixed  step  size  may  reduce  classification 
accuracy; 

(2)  the  storage,  speed,  and*accuracy  are  high¬ 
ly  dependent  on  vector  ordering  and  step 
size  selection. 

Description  of  Method 

The  procedure  to  store  the  library  can  be 
described  as  follows: 

(1)  Arrange  the  N-dimensionat  vectors  in  an 
order  which  minimizes  incremented  dis¬ 
tances  between  the  vectors  (a  possible 
ordering  would  be  by  rows  as  shown  in  Fig. 
5). 

(2)  Code  the  first  N-dimensional  vector 
directly.  Choose  a  step  size  S  for  the 
differential  coding.  Replace  each  subse¬ 
quent  vector  by  its  quantized  difference 
from  the  previous  quantized  vector.  Each 
difference  vector  dimension  is  quantized 
by  rounding  to  an  integer  multiple  of  the 
step  size  S. 

(3)  Replace  each  quantized  N-dimensional  in¬ 
crement  vector  Ly  a  succession  of  elemen¬ 
tary  vectors  (N-dimensional  chain  code) 
with  values  +1,  -1,  or  0  in  each  dimension 
representing  an  increase  of  size  S,  de¬ 
crease  of  size  S,  or  no  change  in  value, 
respect i vely. 

If  aspect  angle  information  is  also  stored, 
the  angle  associated  with  each  incremental 
vector  is  linearly  interpolated  between  li¬ 
brary  vectors. 

As  an  example  consider  the  3-dimensional 

increment  vector  (5, 3, -2)  with  a  step  size  of 

1.  Assume  the  aspect  angle  associated  with 

the  previous  vector  is  8  =  0.5,  8  =  0.5  and 

x  y 

with  the  present  vector  is  8  =  0.0,  e  =  1.0. 

X  7 

The  largest  difference  value  is  5,  therefore 
five  elementary  vectors  will  be  required  to 
represent  this  •  difference  vector.  Table  1 
shows  the  incremental  vector  decomposition. 

The  classification  procedure  is  done  by 
scanning  through  the  library  looking  for  the 
nearest  neighbor  in  N-dimensional  space  to  the 
vector  ••epresenting  the  normalized  unknown  ob¬ 
ject.  A  faster  classification  may  be  obtained 
by  computing  each  distance  by  updating  the 
previous  distance  using  the  incremental  vec¬ 
tors. 

The  step  size  selection  has  a  critical 
effect  on  performance  of  this  storage  pro¬ 
cedure.  As  the  step  size  is  increased  the 
following  effects  are  observed: 

(1)  the  number  of  incremental  vectors  r  .ired 
decreases  and  thus  the  storage  . wquired 
and  classif ication  time  decrease'. 

(2)  the  quantization  error  is  representing  the 
1 ibrary  increase  and  thus  the  classifica¬ 
tion  accuracy  and  aspect  angle  estimation 
accuracy  decreases. 
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Table  1.  Decomposition  into  elementary  vectors 


Step 

Sampled 

Increment 

Vector 

Chain  Code 

Vector 

Angle 

xl 

x2 

x3 

xl 

x2 

x3 

9 

e 

X 

y 

1st 

1 

0.6 

-0.4 

1 

1 

0 

0.4 

0.6 

2nd 

•» 

c. 

1.2 

-0.8 

1 

0 

-1 

0.3 

0.2 

3rd 

3 

1.8 

-1.2 

1 

1 

0 

0.2 

0.8 

4th 

4 

2.4 

-1.6 

1 

0 

-1 

0.1 

0.9 

5th 

5 

3.0 

-2.0 

1 

1 

0 

0.0 

1.0 

Thus  the  step  size  selection  allows  a  trade¬ 
off  between  accuracy  and  the  required  storage 
and  processing  time. 

Library  Vector  Ordering 

Another  important  factor  which  affects 
storage,  speed,  and  accuracy  is  the  order  in 
which  the  library  vectors  are  arranged.  The 
incremental  method  chosen  for  storing  the  li¬ 
brary  can  be  thought  of  as  representing  the 
toroidal  surface  in  N-dimensional  space  with  a 
line  in  N-dimensional  space  which  traces  over 


the  surface.  The  library  vector  ordering 
determines  the  path  that  line  will  take  as  it 
attempts  to  cover  the  surface. 

As  an  example  consider  ordering  the  vec¬ 
tors  in  row  order  as  shown  in  Fig.  5.  The  in¬ 
cremental  vector  locations  when  a  step  size  of 
350  is  used  map  into  the  locations  shown  in 
Fig  6.  Note  that  small  changes  in  9^  angle 

can  cause  large  changes  in  the  N-dimensional 
vector  as  evidenced  by  the  large  number  of  in¬ 
cremental  vectors  at  the  left  center  and  right 
center  of  Fig.  6.  If  column  ordering  were 
used  instead  of  row  ordering,  the  angles  asso- 
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Fig.  5.  Library  views  in  row  order.  Each 
circle  represents  an  aspect  angle  used  in  the 
original  library. 


0X  ANGLE 


Fig.  6.  Library  views  in  row  order  after 
chain  code  representation.  Each  circle 
represents  the  aspect  angle  for  each  chain 
code  increment.  (Step  =  350). 


ciated  with  the  incremental  vector  locations 
are  as  shown  in  Fig.  7. 

One  might  attempt  to  find  the  optimal 
ordering  of  the  library  vectors  so  that  the 
fewest  total  incremental  vectors  are  required. 
The  problem  is  to  find  the  shortest  path 
through  the  Library,  including  each  vector 
once  and  on{.y  once.  This  is  the  wetl  known 
problem  of  the  Traveling  Salesman.  Finding 
the  optimum  solution  is  feasible  only  for 
problems  with  a  small  number  of  nodes.  Howev¬ 
er  many  heuristic  procedures  have  been 
developed  for  approximating  the  optimal  solu- 
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tion.  8.  Golden  et  al.  have  compared  some 
approximate  solution  algorithms  and  published 
an  interesting  review  of  this  problem.  Some 
are  very  simple,  computational  efficient,  and 
give  excellent  approximation. 

The  approximation  we  will  use  is  called 
"nearest  insertion."  The  procedure  is  as  fol¬ 
lows: 

(1)  Consider  all  N-dimensional  vectors  to  be 
nodes. 

(2)  Start  with  a  subgraph  of  node  i  only. 

(3)  Find  node  k  such  that  the  distance  from 


node  i  to  node  k  is  minimal  and  form  sub¬ 
tour  i-k-i. 

(4)  Find  another  node  k  not  in  the  subtree 
closest  to  any  node  in  the  subtour. 

(5)  Find  the  a'rc  (i,j)  in  the  subtour  which 
minimized  d.  .  +  d.  .  -  d.  .  and  insert  k 

1,K  K,J  1,J 

between  i  and  j . 

(6)  Go  to  step  (4)  unless  all  nodes  are  al¬ 
ready  contained  in  the  path. 

It  can  be  shown  that  the  worst  case  path 
length  using  this  procedure  is  less  than  or 
equal  to  twice  the  optimal  length.  The  algo¬ 
rithm  requires  on  the  order  of  N*  computa¬ 
tions.  In  our  implements  each  node  is  select¬ 
ed  as  a  starting  node  for  the  procedure  and 
the  best  result  retained.  Therefore,  the  com¬ 
putational  complexity  becomes  D(N^). 

8.  Results 

As  a  test  of  the  methods  proposed  in  Sec¬ 
tion  7  an  experiment  was  run  using  6  airplane 
libraries  (mirage,  mig,  phantom,  B-57,  F104 

and  F105).  These  are  the  same  shapes  used  i 
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Fig.  7.  Library  views  in  column  order  after 
chain  code  representation.  Each  circle 
represents  the  aspect  angle  for  each  chain 
code  increment  (Step  =  350). 


References  5,  10,  and  13.  For  all  classifica¬ 
tion  experiments,  we  considered  an  unknown  set 
composed  of  600  views  (each  are  a  128  x  128 
grid),  100  for  each  type  of  airplane.  The 
normalized  Fourier  descriptor  was  calculated 
and  compared  to  the  library  entries  which  were 
stored  in  the  N-dimensional  chain  code  format 
described  in  Section  6. 

Shown  in  Table  2  are  some  experimental 
results.  For  example  in  entry  (1)  of  Table  2 
row  ordering  of  each  original  library  (see 
Fig.  5)  was  used  prior  to  converting  to  a 
12-dimensional  code.  The  unknown  airplanes 
were  classified  correctly  with  an  accuracy  of 
91X.  The  median  estimated  angle  error  was 
.075  radians  in  0^  (along  the  rows)  and  .082 

radians  in  0^  (along  the  columns).  The 

storage  required  for  the  12-dimensional  chain 
codes  for  the  6  libraries  was  43.0  k-bytes. 

As  the  step  size  of  the  chain  code  is  in¬ 
creased,  the  quantization  error  in  the  library 
vectors  increases,  but  the  storage  required 
decreases.  Thus  the  classification  accuracy 
is  decreased  but  the  storage  requirements  and 
the  processing  time  decreases.  This  is  shown 


Fig.  8.  Representative  contours  of  the  six 
airplanes  used  for  the  experimental  results. 
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Table  2  Experimental  Results. 

Unless  otherwise  stated:  row  ordering,  absolute  distance  criterion,  N  *  12  dimensions,  each  ori¬ 
ginal  library  has  143  views. 


Experimental 

Test 

(S  =  step  size) 

Classification 
Accuracy  (X) 

Median  Angle 
Error  (radians) 
e  e 

_ * _ L 

Storage 

Required 

(k-bytes) 

(1) 

S  =  100 

91.0 

Hilim 

43.0 

(2) 

S  =  200 

89.8 

.080 

.082 

21.9 

(3) 

S  =  300 

88.8 

.095 

.082 

14.7 

(4) 

S  =  400 

84.3 

.104 

.083 

11.3 

(5) 

S  =  500 

81.5 

.109 

.083 

9.1 

(6) 

column  ordering,  S  =  200 

89.0 

.131 

.048 

28.9 

(7) 

column  ordering,  S  =  500 

81.2 

.144 

.066 

12.0 

(8) 

optimal  ordering,  S  =  200 

89.0 

.133 

.078 

14.4 

(9) 

optimal  ordering,  S  =  500 

82.7 

.144 

.084 

6.4 

(10) 

mse  distance,  S  =  200 

86.8 

.071 

.082 

21.9 

(11) 

augmented  lib  (525  views). 

S  =  200  95.8 

.051 

.046 

49.4 

(12) 

aug.  lib,  N  =  16,  S  =  200 

97.3 

.043 

.046 

60.0 

in  entries  (2)  -  (5)  of  Table  2. 

If  column  ordering  of  the  original  library 
is  used  (see  Fig.  6)  the  accuracy  and  storage 
requirements  both  degrade  compared  to  row  ord¬ 
ering  (see  entries  (6)  and  (7)  in  Table  2). 
Thus  it  is  concluded  that  row  ordering  results 
in  more  average  correlation  of  the  sequential 
12-dimensional  vectors  than  does  column  order¬ 
ing. 

Shown  in  entries  (8)  and  (9)  of  Table  2  are 
the  results  of  using  the  "nearest  insertion'* 
vector  ordering  described  at  the  end  of  Sec¬ 
tion  7.  Although  the  storage  requirements  are 
better,  the  accuracy  decreases  somewhat.  Upon 
close  inspection  of  the  details  of  this  order¬ 
ing,  we  have  found  that  the  views  of  the  air¬ 
plane  for  quite  different  angles  may  be  very 
similar.  The  optimal  ordering  thus  tends  to 
have  large  jumps  in  aspect  angle  which  results 
in  poor  interpolation  and  degrades  classifica¬ 
tion.  Improvements  are  now  underway  to  limit 
the  optimal  ordering  to  selected  regions  of 
aspect  angles  and  we  feel  classification  and 
angle  accuracies  will  improve  significantly. 

Entry  (10)  shows  the  results  of  using  a 
mean-square  error  distance  criterion  instead 
of  mean  absolute  distance.  Accuracy  is  some¬ 
what  degraded.  Me  feel  the  mean-square  cri¬ 
terion  tends  to  weight  the  targe  power  low 
harmonic  terms  too  much  while  the  airplane 
type  is  discriminated  using  the  lower  power 
higher  harmonic  terms. 

If  the  original  library  is  increased  to 
25  x  21  instead  of  13  x  11,  the  storage  re¬ 
quired  increases  as  shown  in  entry  (11)  but 
the  classification  accuracy  increases  dramati¬ 
cally.  Further  increase  in  classification  ac¬ 
curacy  can  be  obtained  by  retaining  16  numbers 


for  each  library  entry  instead  of  12.  This  is 
shown  in  entry  (12)  of  Table  2. 
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